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Abstract. In this paper we study Tikhonov regularization for the stable so- 
lution of an ill-posed non-linear operator equation. The operator we consider, 
which is related to an active contour model for image segmentation, is contin- 
uous, compact, but nowhere differentiable. Nevertheless we are able to derive 
convergence rates under different smoothness assumptions on the true solu- 
tion by employing the method of variational or source inequalities. With this 
approach, we can prove up to linear convergence with respect to the norm. 
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1. Introduction 

The study of convergence rates within the context of solutions of ill-posed prob- 
lems is concerned with the question of how good one can approximate the true 
solution of an equation given some noisy data of known noise level 6. Assume that 
A : X — > Y is a bounded linear operator between two Hilbert spaces X and Y and 
we are trying to solve the linear operator equation Ax — y s with noisy right hand 
side y 5 satisfying \\y — y s \\ < 6 for some "true" data y 6 Y. Then the classical 
Tikhonov approach consists in the minimization of the functional 

(1) \\Ax-y s fy + a\\xf x 

with some regularization parameter a > 0. If 6 and a tend to zero in such a way 
that S 2 /a — > 0, then the minimizers x s a of ([T]) will converge to the "true" solution 

x f = A f y 

with A< being the Moore-Penrose inverse of A [11]. Moreover, it is possible to 
estimate the quality of the approximation, that is, the difference — a^H , provided 
one has additional a-priori knowledge about smoothness properties of x': If x* = 
(A*A) w ijj for some ui € X and < v < 1, then one can show for the parameter 
choice a ~ that we have (TTJ Corollary 3.1.4] 

(2) 11^-^11=0(5^) as <S->0. 

In the case of a non-linear operator equation F(x) — y, the classical approach 
to the derivation of convergence rates uses a linearization of F at the true solution 
x^. Then, analogous convergence rates as in the linear case are derived under the 
assumption of a source condition of the form 

(3) (F'(x^)*F'(x^)) v u; — xK 
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In the proofs of these convergence rates, however, it is necessary to introduce ad- 
ditional assumptions coupling the non-linear functional F and its linearization at 
it. 

An alternative approach to non-linear problems has been suggested in [16 . In- 
stead of linearizing the problem and using source equalities, they proposed to use 
a single source inequality of the form 

(4) ||i - x^f x < \\x\\\ - ||it||| + h\\F(x) - F{xt)\\ Y . 

With this assumption, they were able to prove a convergence rate in the norm of 
order 0(\/~5). This corresponds to the case v — 1/2 in ([2]); indeed, in the linear 
case, the source inequality (j4|) is equivalent to the source condition it = (A* A)zw 
(see [I!]). Moreover, the approach has been generalized in [2j [7] to intermediate 
convergence rates, roughly corresponding to the source condition it — (A* Ayuj 
with < v < 1/2. 

Apart from the possibility of treating non-linear operators, source inequalities 
have the advantage that they do not require a Hilbert space structure of the set 
X; indeed, the results in [16] were formulated in a Banach space setting ([7] does 
not even require linear spaces). One problem in Banach spaces is that fractional 
powers of a linear operator are not defined any more; a source condition of the 
form ((3]) cannot even be formulated. There is, however, a different approach to 
convergence rates of general order for non-linear problems in Banach spaces by 
means of approximate source conditions (see |12 | I13 [ H4 l 115]). Still this approach 
requires a prior linearization of the operator. 

While theoretically the approach by source inequalities can be used for deriving 
convergence rates in general non-smooth situations, most cases where it is been 
applied were differentiable settings with a weak coupling between the operator and 
its linearization. In this paper we will show by means of a concrete example that the 
method of source inequalities can also be used for deriving convergence rates under 
meaningful conditions for a non-diffcrcntiable operator. We consider a situation 
that is closely related to an active contour model for image segmentation, where 
the operator, in its natural Hilbert space setting, is Holder continuous but nowhere 
differentiable. Nevertheless it is possible to derive convergence rates up to order 
0(S) under the assumption of sufficient smoothness of the true solution (in the 
sense that it lies in a certain Sobolev space). 

2. Setting 

Denote by H^_(S 1 ) the non- negative cone in the Sobolev space H 1 (S 1 ). That is, 
7 £ HUS 1 ) if and only if 7 e H^S 1 ) and j(t) > for every t £ S 1 . If 7 £ H^S 1 ), 
we denote by 

hyp(7) := {(t,x) £ S 1 x K> : < x < 7 (i)} 
the hypograph of 7. That is, hyp(7) consists of all points lying between the i-axis 
and the graph of 7. We now define the mapping F: H^S 1 ) — > L 2 (S' 1 x M>o), 

F(u) := Xhy P ( 7 ) 

with 

'0 if7(t)>i, 
1 if < j(t) < x, 
which maps the curve 7 to the characteristic function of its hypograph. 



Xhyp(7) (j't l) 
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Note that every function 7 <E ^[^(S 1 ) is bounded and continuous, and therefore 
hyp(7) C S 1 x K>o is a (non-empty) compact set. In particular this implies that the 
mapping F indeed takes values in the space L 2 (S' 1 x R>o). Moreover it is obvious 
that F is injective. 

We now consider the problem of solving the equation ^(7) = u for some given 
u G Ran(_F) in a stable way. In addition, we assume that, instead of u, we are only 
given noisy data u 5 € L 2 (S' 1 x R>o) satisfying \\u s — u\\l2 < 5, where S > is some 
known noise level. Here we do not assume that the noisy data lie in the range of F 
or even that they are the characteristic function of some (Borel) set. 

In order to find an approximate solution of the equation 

(5) F( 7 ) = u 5 , 

we apply classical (quadratic) Tikhonov regularization using the squared homoge- 
neous first order Sobolev norm as regularization term. That is, we consider the 
Tikhonov functional 

r«(7;«'):= 11^(7) -«*||i a +a||7lli> 
and define, for a > 0, an approximate solution of ([5} by 

7* :e argmin{r Q ( 7 ;« 5 ) : 7 € ^K^ 1 )}. 

Relation to Active Contours. One of the basic problems in image processing is 
segmentation: Given a possibly noisy image, interpreted as a function u € L 2 (R 2 ), 
the task in its simplest form is to partition this image into two regions, one repre- 
senting objects in the image, the other representing the background. To that end, 
Chan and Vese [3] proposed a variational method based on the assumption that the 
objects (which are assumed to form a simply connected region) share a common 
grey value c\ , and also the background has a uniform grey value C2 ^ c\ . Without 
loss of generality we may assume that c-i — 1 and c\ = 0. Then the method pro- 
posed in [3] consists in minimizing, for some regularization parameters ai, ai > 0, 
the functional 

(6) 1(7) := ai length(7) -(- a2 area(7) + / (1 — u(x)) 2 dx + / u(x) 2 dx. 

•Ant (7) ^K 2 \int(7) 

Here length^) denotes the length of the curve 7, int(7) C M 2 denotes the set 
of all points lying inside the curve 7, and area(7) the size of the area inside the 
curve. This functional is minimized over the space Emb(5 1 ; R 2 ) of all closed curves 
continuously embedded in R 2 . 

Defining the function F: Emb(S' 1 ;R 2 ) -> L 2 (R 2 ), 

F(7) ■= Xint( 7 ), 

one can write the functional cquivalently as 

1(7) = ai length(7) + a 2 area(7) + ||F( 7 ) - u||| 2 . 
Moreover, for 7 sufficiently smooth, one can write 

length( 7 )=/ m)\dt, area(7) = i/ , ( 7 (i), 7 (t)) di. 

Thus the Chan- Vese model ([6]) can be brought in a quite similar form as the 
problem we consider in this paper, which, conversely, can be interpreted as a sim- 
plified active contour model for images on a half-cylinder. A major reason for our 
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simplification is the fact that it allows us to work in a Hilbert space setting. In 
contrast, the "natural" setting for the Chan-Vese model would rather be the Ba- 
nach manifold of Lipschitz embeddings of S 1 in R 2 modulo reparameterizations of 
S 1 (note that ^(7) is invariant with respect to reparameterizations of 7; in order 
to obtain any uniqueness result it is therefore necessary to factor out reparameter- 
izations). 

3. Well-posedness of the Regularization Method 

We now prove that quadratic Tikhonov regularization yields a well-posed regu- 
larization method for the functional F. To that end we mainly have to investigate 
the continuity properties of F. 

Lemma 3.1. The mapping F is Holder continuous of degree 1/2 as a mapping 
from L 2 ^(S 1 ) to L 2 (S 1 x R>o). In particular, it is weakly continuous and compact 
as a mapping from H^(S 1 ) to L 2 {S 1 x M>o)- 

Proof. Let 71, 72 G L+tS" 1 )- Then 

||F(7i) -F(<y 2 )\\h> = / / IXhy P ( 7l )(i,z) - Xhyp( l2 )(t,x)\ 2 dxdt 

JS 1 JR>0 

|7i(*) - 72(01 dt = Il7i - 72k 1 < V^7r||7i - 72IU 2 , 

which proves the Holder continuity of F as a mapping from L 2 i _(S 1 ). The weak 
continuity and compactness of F as a mapping from H^S 1 ) now follow from the 
fact that i? 1 (S' 1 ) is compactly embedded in L 2 (S' 1 ). □ 

Theorem 3.2. The following hold: 

(1) For every u e L 2 {S 1 x K>o) and a > the functional T a (-;u) attains its 
minimum in H^S 1 ). 

(2) Assume that u G L 2 (S 1 x K> ) and a > 0. Let {-u (fc) } fce N C L 2 (S l x M> ) 
be any sequence converging to u, and let 

7 (fe) G argmin{ra(7;w ( ' £) ) : 7 G ^(S 1 )}. 

Then there exists a sub-sequence {T^'ljeN converging strongly in H 1 (S 1 ) 
to some 

7 G argmin{r a (7;M) : 7 G HKS 1 )}. 

(3) Assume that u — Ffa') for some j> G H^(S 1 ). Let {5k}keN C K>o and 
{afcjfeGN C M>o be sequences satisfying 6k — > 0, ctk — > 0, and S 2 /ak — > 0, 
ana letu 8 " G L 2 (S 1 x M>o), k G N, satisfy \\uk — «||l 2 < $k- Let 

7 (fc) := li\ G argmin{T Qfc (7;^) : 7 G ^(S 1 )}. 

T7ien the sequence {"f^^^&i converges strongly in H 1 (S 1 ) to 7^. 

Proof. This follows from standard results in the theory of (non-linear) Tikhonov 
regularization, which can be found, for instance, in [H 119] . The only (ever so 
slight) complication is the fact that the regularization term uses the homogeneous 
Sobolev norm, which is not coercive on H 1 ^ 1 ). The coercivity of the functional 
7~ a , however, follows immediately from the chain of inequalities 

||F( 7 ) - uf L2 < 2||F( 7 )|| i2 + 2||u||| a = 2|| 7 || £ x + 2|M| 2 L2 < 2|| 7 |U> + 2||u||£ a . 
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□ 

Remark 1. Because F is injective, it follows that the equation F(-f) — u can have 
at most one solution. Still it is possible that the optimisation problem Ta(^;u s ) — > 
min can have multiple solutions, even for 6 > arbitrarily small. 

Consider for instance the situation where 7^(t) = 1 and F(^) — Xs 1 x[o,i]' 
Define for 5 > with S 2 < tt the function 

{1 if < x < 1 - <S 2 /tt, 

~ if \x-l\ <5 2 /ir, 

if x > 1 + 5 2 /ir. 

Then \\u s ~u\\ L 2 = S. 

Now let 7 e H^S 1 ) and define 

J+{5) := {teS 1 : 7 (t) > 1 + 5 2 /tt}, J- (5) := {t e S 1 : 7 (t) < 1 - S 2 /ir}. 
Then it is easy to see that 

11^(7) " u s \\ 2 L 2 =5 2 + [ 1 (t)-l-5 2 /irdt+ [ l- 7 (<)- ( 5 2 /7Tdt. 

Jj+(8) Jj-(S) 

Therefore ||-F(7) — u '|lia 1S minimal (with value S 2 ) if and only if \j(t) — 1| < 5 2 /n 
for every t. Thus it follows that 

T a (r,u s ) >5 2 

for every 7 g H+tS 1 ), and equality holds if and only if 7 is a constant function 
of the form 7(<) = 1 + c for some |c| < S 2 /tt. In other words, every such function 
is a minimizer of the Tikhonov functional with data u s and any regularization 
parameter a > 0. 

In particular, this example shows that the stability result in Theorem 13.21 (item 
(2)) really requires the formulation in terms of sub-sequences; it can happen that the 
sequence {7^ fe ^}fc e N itself does not converge. Note, however, that it is also possible 
to formulate the stability result in terms of set convergence. Such an approach has 
for instance been used in [ID] . 

4. Convergence Rates 

In this main section of the paper we will discuss the derivation of convergence 
rates, that is, quantitative estimates for the difference between the regularized 
solution 7* and the true solution 7^ in dependence of a and 8. 

Throughout this whole section we assume that the equation ^(7) = u has a (nec- 
essarily unique) solution 7^ € H^(S 1 ). Moreover we denote by 7* any minimizer 
of the regularization functional 7^(-;t*) for any noisy data u s e L 2 ^ 1 x M> ) 
satisfying \\u s — u\\ L 2 < S. 

Classically, convergence rates for quadratic Tikhonov regularization of non-linear 
operators on Hilbert spaces have been derived for sufficiently smooth operators un- 
der the assumption that the true solution 7^ is contained in the range of the adjoint 
of the derivative of the operator F at 7^ . In addition it is necessary to compensate 
for the non-linearity of the operator F by assuming additional regularity properties 
of its derivative F' . One example is the following theorem taken from [15] . Similar 
results can also be found in [4] [5] . 
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Theorem 4.1. Assume that F is Gateaux differentiable in a neighborhood of ^ 
and that there exist lu £ L 2 (S' 1 x R>o) and c > with c||cj||x,2 < 1 such that 

and 



\\F(l) ~ Ftf) - F'( 7 t)( 7 - 7 t)|| L2 < £|| 7 - 7 t||2 ffli 

for every 7 £ H^S 1 ) sufficiently close to 7 ^. Then we have for a parameter choice 
a ~ S the rate _ 

\\7 S a-7 f \\m=0(VS). 

The crucial assumption in this theorem is the source condition 7^ £ F'(j^)*uj, 
which links the non-linear problem to the much better understood linear theory. 
The additional assumptions mainly guarantee that the gap between the non-linear 
and the linearized problem is sufficiently small. Still, the whole approach for the 
derivation of convergence rates relies on a linearization of the operator F, which in 
our situation is not easily possible, as the next result shows. 

Lemma 4.2. The mapping F: H\(S ) — > L 2 ^ 1 x R>o) is nowhere differentiable 
in the interior of its domain. 

Proof. First note that the interior of H^(S 1 ) consists of all functions 7 £ iJ 1 (S' 1 ) 
with j(t) > for every t £ S 1 . 

Let therefore 7 £ H]_(S 1 ) satisfy j(t) > everywhere, and let a £ H 1 (S 1 ). 
Then, as shown in the proof of Lemma 13 . 1 1 we have 



\F(j+sa)-F(j)\\l 2 = f \sa(t)\dt 



Now assume that a ^ 0. Then there exist a non-empty open subset J C S 1 and 
e > such that |er(i)| > e for every t £ J. Consequently 



|F( 7 + S a)-F( 7 )f f2 = | S | / \a(t)\dt>e\s\\J\. 

s 1 



Thus 



lim -||F( 7 + sa) - F( 7 )|| r2 > lim \ = +00, 
s->o+ s UJj s^o+ V s 



which proves that F does not even have a one-sided directional derivative in any 
non-trivial direction. □ 

For the next result — the main theorem of this paper — recall the definition of the 
fractional order Sobolev spaces W s ' q , s £ K, 1 < q < 00, and the homogeneous 
Sobolev spaces W^ q , which can for instance be found in pQ. 

Theorem 4.3. Assume that 7+ £ W s ' q (S 1 ) with 1 < s < 2 and 2 < q < +00. 
Then there exist constants c\, C2, C3 > such that 

cib-yil^ 

< hWll - iMll^ + c 2 ||F( 7 ) - Ftf)\fe-° + c 3 ||F( 7 ) - F( 7 t)£^i 

for every 7 £ H+IS 1 ). 

J/ 7 t g W'-^iS 1 ) with 1 <s<2, then 



(7) Cl ||7 - 7 + ll^ < IMIij - \tff Hi + C2\\F(j) - F( 7 t)|| 



L 2 
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Proof. We only show the assertion for the more complicated case q < oo. 
First note that 

Il7-7 t l| 2 = ll7l| 2 -||7 t l| 2 -2(7 t ,7-7 t >- 

In order to estimate the last term on the right hand side of this equation, we will 
use the abbreviation 

er := 7 — 7^. 
Because 7* € W+^S 1 ), it follows that 

\(¥,o-)\<hH\wr\\<r\\w, 



2-3,9, 



Moreover the interpolation inequality for Sobolev functions implies that there 
exists some C\ > only depending on s and such that 

11,-rll << C II ^-11 S_ 1 II ^-11 2 — S ^ SI II II s — 1 II _||2 — s 

\\a\\ w 2-s, q , < Ci||o-|| L „„ ||0"|| i,„ < Oi||o-|| i9 , \\<T\\ H i ■ 
The last inequality holds, because q > 2 and therefore < 2. Now note that 

IMU«. < IkllfelkllL 

< Ikll^dHUi + lkll^O" 
<IHUi + |H|§H^ ol . 



Setting C2 := 2C\ ||t II w^' a ; we thus obtain the estimate 

2|<7 + ,*>| <c , 2|H|'T 1 IH|^r + c' 2 |H|§ 1 |H|^ 1 



Now we use Young's inequality ab < ^a p + -^-b p » , which implies with p = 2/s and 
p* = 2/(2 - s) that 

CaH^IM&r ^ f ^1 IMIlf + (1 - §)lMI*j 

and with p = and p* = ^rj^r 



Thus we obtain with the constants 

(2g-l)(*-l) 



Cl = 



2? 



_ ga - 8 + l 'j^ 
C3 - 2g ^ : 



the estimate 



n ||- - Vll 2 ^ < IMI^ - Il7 f ll^ + cah - + c 3 ||7 - 7 f lll? 



Now the assertion follows from the fact that ||F(7) - F(rf)\\ 2 L2 = ||7 — 7^||l 1 - ^ 
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Corollary 4.4. Let the conditions of Theorem \4-3\ be satisfied and let 

$(t) : = C2t 2£ ^ + oste- 
in case q < oo and 

$(t) := c 2 t 2 ^ 

in case q = oo. 

(1) If q < oo and s < 2 let ^! be the convex conjugate of the strictly convex 
mapping t i— > (2i) . TTien 

Cl ||7^-7 t ll^<^ + ^ 2 ) + ^. 

(2) If q = oo and s = 2 we have 

cihi-lH 2 Hl <^ + 2c 2 6 2 
whenever a < l/(2c2). 
Proof. This is a direct consequence of Theorem l4.3l and [71 Theorem 3.1]. □ 

Corollary 4.5. Lei i/ie conditions of Theorem \4-3\ be satisfied. 

(1) If q < oo or s < 2 we have for the parameter choice a ~ S i* the 
convergence rate 

ll7£-7 + lk =0(5&). 

(2) If q = oo and s = 2 we have for a constant parameter choice a < \j(2c2) 
the convergence rate 

\\li-l X \\Hl=0{H). 

Proof. See [7, Corollary 3.1]. □ 



Remark 2. Note that the claim of Theorem 14.31 in the case s = 2 and q — oo can 
be shown much more easily with elementary computations. Indeed, in this case one 
has 



IIT-^IIL = II7IIL " ll7 f HLx -2 / 7 t (*)(7(*)-7 t (*)).d* 

ooo j sl 

= ll7l| 2 fff; -||7 t ll^+2 / 7 t W(7(*)-7 t (*))* 
Js 1 

<ll7l| 2 ff( ;-||7 t ll^+2||7 t l|L~||7-7 t IU 1 

= ll7l| 2 ff i-||7 t ll^+2||7 t ||L»||F( 7 )-F( 7 t)|| 

Thus (UJ holds with c x = 1 and c 2 = 2\\^\\ L oa. 



2 

L 2 ' 



5. Further Aspects 

5.1. Differentiability. In Lemma T4. 2 1 we have shown that F is nowhere differen- 
tiable when regarded as a mapping from H^_(S 1 ) to L 2 (S 1 x R>o). We now discuss 
possible different settings, where F has better regularity properties. 

Because the range of F consists only of characteristic functions with compact 
support, it follows that F can be regarded as a mapping into any space L P (S 1 x K>o) 
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with 1 < p < oo. Moreover the same computation as in the proof of Lemma 13.11 
shows that for every 1 < p < oo we have 



Thus, seen as a mapping F: i 2 (5 1 ) — > L p (S' 1 ;R>o), the mapping F is Holder 
continuous of degree 1/p; in case p = 1, this shows that F is Lipschitz continuous. 
Note that for p = oo, the mapping F is discontinuous everywhere. Concerning the 
differentiability of F, however, this change of the target space does not really matter. 
The same argumentation as in Lemma T4.2I shows that F is nowhere differentiable 
as a mapping into L P (S' 1 x R>o) with 1 < p < oo. 

For p — 1 the situation is slightly different: If 7 6 Hi_(P) satisfies j(t) > for ev- 
ery t and a £ H 1 ^) is fixed, then the family of functions {^(F(j + s<r)—F(-f))} s>0 
forms abounded subset of xR> ). Still, the limit lim s ^. + |(^(7+ scr ) — 

does not exist in x R>o) and thus, again, F is nowhere directionally differ- 

entiable. 

Finally, one can regard F as a mapping from H^_(S 1 ) to the space MiS 1 x M>o) 
of finite Radon measures on S 1 x M>o. Because x R>o) is isometrically 

embedded into A^S 11 x K>o), this does not change the Lipschitz continuity of F. 
However, the difference quotients -(^(7 + sa) — F(~f)) of F now have a limit in 
^(S 1 x R>) with respect to the weak* topology: One has 



where 7^(s£ 1 ) denotes the push forward of the one-dimensional Lebesgue measure 
C 1 weighted by the functions a via the mapping 7. That is, 



for every h £ C^S 1 x R> ). 

The considerations above show that the setting of the problem can be modified 
in such a way that the operator F becomes differentiable. Thus it might still be 
possible to derive convergence rates in a more classical way via a linearization of 
F. There are, however, several difficulties: First, the extension to ^(S 1 x M>o) 
leaves the Hilbert space setting in favor of a more complicate Banach space setting 
(with non- reflexive spaces). In such a setting, non-standard convergence rates have 
been derived via linearization ideas in |14j by the technique of approximate source 
conditions. Second, even the Banach space setting might not be the correct one, 
as the operator F is only differentiable with respect to the weak* topology and not 
the norm topology. Thus it might even be necessary to derive results in general 
locally convex spaces. Finally note that after changing the topological structure, 
the similarity term is not the squared norm on the target space any more. Thus it 
might even be necessary to turn to results for convergence rates with more general 
similarity terms (see for instance [6l El [18] ) . 

5.2. Differentiability of the Regularization Term. Interestingly, the regular- 
ization functional itself has potentially much better smoothness properties than F. 
In order to see this, define for u £ i 2 (5' 1 x R>o) the functional S: H^S 1 ) — > M>o, 



\\ F (li) - F (l2)\\ P L p = II71-72IU1 < V2n\h - 12\\l 2 - 



(F( 7 + sa) - Ffr)) ^* 7 # (^ 1 ) 




Suh) := \\F( 7 )-uf L2 . 
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Then it is easy to sec that S u is Frechet diffcrcntiablc whenever u £ L 2 (S' 1 x R>o)n 
C(S 1 x M>o) with Frechet derivative 



L Js 1 J 
In addition, for u — F(-y) the mapping S, while not differentiable, possesses one- 
sided directional derivatives of the form 



The specific form of this one-sided derivative might possibly yield an explana- 
tion of the linear convergence rate that has been shown in Corollary [13] for the case 
where 7^ £ W 2 ' 00 (S 1 ), as such a rate is not possible for smooth operators using the 
squared Hilbert space norm as a regularization term. For quadratic regularization 
of Frechet differentiable operators, it has been shown that the best possible con- 
vergence rate in non-trivial cases is of order 0(<5 2 / 3 ) (see [17]). Linear convergence 
rates, however, have been derived recently for Tikhonov regularization with non- 
smooth regularization terms. The first result in this direction was [9], where linear 
rates have been derived for regularization on the sequence space £ 2 with the ^-norm 
as regularization term under the assumptions of sparsity of the true solution and a 
certain source condition. This result has also been extended in [5] to more general 
positively homogeneous regularization terms. The basis of all these results is the 
non-smoothness of the regularization term; the rates derived in this paper indicate 
that a sufficient non-smoothness of the operator F to be inverted can have similar 
effects. 



We have shown in this paper that the approach of source inequalities for the 
derivation of convergence rates for Tikhonov regularization can be applied to non- 
linear problems where approaches based on linearization are bound to fail. Our 
paradigm was a functional related to the Chan-Vese active contour model for image 
segmentation, which was shown to be continuous but nowhere differentiable in its 
domain. Still, the approach by source inequalities allowed us to derive convergence 
rates under reasonable (and easily interpretable) smoothness assumptions on the 
true solution. Surprisingly, the convergence rates do not obey the classical bound 
of 0(<5 2 / 3 ), which is known to be best possible rate for quadratic regularization of 
differentiable functionals. Instead, for a sufficiently smooth true solution, we were 
able to obtain a rate of order 0(S). One possible explanation for this exceedingly 
good behaviour is a connection to sparse (or i 1 ) regularization: For noise free 
data the similarity term in our problem has precisely the same behaviour as the 
^-regularization term. 
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